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1. Introduction 



Recent calculations in both the meson [||, ^ and baryon [^] sectors within the twisted mass 
formulation of lattice QCD [Q] have yielded accurate results on a wide range of observables of 
immediate relevance to experiment and phenomenology This is a consequence of the fact 
that the simulations of the European Twisted Mass Collaboration (ETMC) cover a well chosen 
range of parameters, allowing controlled continuum and chiral extrapolations and thus making a 
reliable connection to the physical regime possible. In the present work we apply the twisted mass 
framework to the calculation of nucleon form factors. 

Electromagnetic and axial form factors of the proton and the neutron are fundamental quan- 
tities that yield information on their internal structure such as their size, magnetization and axial 
charge. They have been studied experimentally for a long time with steadily increasing precision, 
the latest generation of experiments reaching relative precisions of down to about 1% The 
first lattice calculations using dynamical fermions appeared only recently 10, ^ and the available 
pion masses were rather large, often above 400 MeV. 

Our goal is to calculate matrix elements of the form {N{pf,Sf)\i&^\N{pi,Si)), where \N{p,s)) 
is the nucleon ground state, with the nucleon having momentum p and spin s. is either the 
electromagnetic or the axial current. Using, in addition, the pseudoscalar current will enable us to 
check phenomenological consequences of chiral symmetry such as the Goldberger-Treiman rela- 



tion [10|. The electromagnetic matrix element can be expressed in terms of two Lorentz invariant 
form factors that depend on the momentum transfer squared only. In Euclidean space-time the 
decomposition is 



{N{pf,Sf)\V^mN{pi,Si)) =u{pf,Sf) 



u{pi,Si 



(1.1) 



where Q = Pf — pt, CJ^v = '[7ai) 7v]/2, m is the proton mass and u{p,s) is a solution to the free Dirac 
equation with mass m. Instead of the Dirac and Pauli form factors Fi and F2 the matrix element 
can be expressed in terms of the electric and magnetic Sachs form factors 

l2 



4m' 



Gm = F,{Q^)+F2{Q^) 



(1.2) 



Similarly the axial current matrix element can be written in terms of the form factors Ga and G„, 



{N{pf,Sf)\A^{0)\N{pi,Si))=u{pf,Sf) 



u{pi,Si). 



(1.3) 



2. Calculation details 



2.1 Wilson twisted mass QCD 

For our calculation we use a Wilson twisted mass fermion action (tmQCD) at maximal twist 
angle together with a tree-level Symanzik improved gauge action. More precisely, the fermion 
action is given by 



1- 



(2.1) 
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with lattice spacing a, the covariant forward and backward lattice derivatives V and V* and mass 
parameters nicru and pL. The Pauli matrix T3 acts on a doublet = {u,d) of two light quarks. 
Automatic 0{a) improvement occurs at maximal twist [11|, which here is realized by tuning the 
value of merit so that the PCAC mass vanishes. The gauge action is 

^^ = |lf| I (l-ReTr[[/i;,V])-i^ I (l-ReTr[[/i,;;2,])V (2.2) 

where U^^^y denotes the plaquette term, U^^\, the rectangular 1x2 Wilson loops and j8 is the 
inverse bare coupling. Further details concerning the simulation and in particular the tuning to 
maximal twist can be found in Ref. 1 12]. 

2.2 Correlation functions 

On the lattice the form factors are extracted from dimensionless ratios of correlation functions. 
We measure the two- and three- point functions 



G{q,t) = £ e-'^/-5r°„ {Utf,Xf)Jp{0)) , (2.3) 
G^,{^\q,t) = £ e"?r^„ {Ja{tf,Xf)ir^,it,x)Jpm , (2.4) 

where = 1 /4[1 + 70] if v = and = 1 /4[1 + JohsYv if V = 1,2 or 3. and / is the proton 
interpolating field. In tmQCD at maximal twist the standard interpolating field reads 

J{x) = -^[1 + /rsle"*" \u'''^{x)'^r5d^{x)] u'\x) , (2.5) 
V2 L J 

with denoting the charge conjugation matrix. To enhance the overlap of J with the proton ground 
state, the quark fields entering Eq. (|2.5[) are smeared. 



u"{t,x) =Y^F'"\x,y-Uit))u\t,y), 
y 

F = {1 + aHf, H{x,y;U{t)) = £ (Ui{x)5,,y^j + uj{x-i)5,^y+, 



(=1 

In addition, the spatial gauge links entering the hopping matrix H are APE-smeared. Good choices 
for the parameters a and N were determined in Ref. [^. 

In this work, we restrict ourselves to the axial- vector and vector currents for the operators ff^ 
inserted into the three point functions. We use the local currents 

(x) = uj^ 75M - dy^ y^d (2.6) 
{x) = uy^u- dy^ ysd, (2.7) 

as well as the symmetrized, conserved vector current 

<W = + where (2.8) 

jUx) = q{x + fL)^[y^,+l]U^'{x)q{x) + q{x)^[y^ - l]U^{x)q{x + jl) . (2.9) 
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While the two-point functions can be calculated using a forward propagator from a fixed 
source, the evaluation of three-point functions is more involved. In correlators containing the 
isovector operators, disconnected diagrams are zero up to lattice artifacts, and can be safely ne- 
glected as we approach the continuum limit. The connected diagrams are calculated using sequen- 
tial inversions through the sink []T3|]. This means that at a fixed source-sink separation we are able 
to obtain results for all possible momentum transfers and insertion times as well as for any operator 
0'^ , with a two sequential inversions per choice of the sink. In this work we use two different sinks, 
one optimized for the electric and one for the magnetic form factor The latter is also suitable 
for the axial form factors. 



2.3 Ratios 

In ratios of three- and two- point functions normalization factors and the leading exponential 
dependencies on the insertion time cancel. From fits to the resulting plateaus, the form factors can 
be extracted. Different ratios can be considered and we here give two examples: 



G{0,2{tf-t))Giq,2it-ti)) 



or 



(r, q,t) / G{q,tf -t)G{0,t)G{Q,tf) 



(2.10) 



(2.11) 



G{0,tf) \ G{0,tf-t)G{q,t)G{q,tf) 
As shown in Fig. [TJ, both lead to compatible plateaus but to different statistical errors. We use 
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Figure 1: Left panel: Comparison of the ratios given in Eqs. (2.10) (fiUed symbols) and (2.11) (open 
symbols) for source type given in Eq. (2.13) and a few representative momentum combinations. Right 
panel: Electric Sachs form factor extracted with a sink source separation of 12fl and 14a. 



the ratio of Eq. (2.11) for the final analysis, which turns out to be superior because it does not 
contain potentially noisy two point functions at large separations and because correlations between 
its different factors reduce the statistical noise. This being most evident at zero momentum with 
the conserved current, where noise cancels completely. 
Once the plateau values 



n(r,^) = lim lim R{r,q,t) 



(2.12) 
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are estimated, the form factors are obtained from the following combinations 

n^(r°,^) = + + Ge{Q^) (2.13) 

k 

U\T\q)+U\T\q)+U\T\q) = ^Y^ej.iqjdij Gm{Q^) (2.14) 



jkl 

n"{r\q)+n'\r\q)+n''{r\q) = ^[{q,+q, + q,)lLGp{Q^)-{E + m)GA{Q')].(2.\5) 



E is the energy of a proton with momentum pj, m its mass and the constant ^ = y E{E+m) 
from the normalization of lattice hadron states that we use. Note that the sequential source needed 
in the calculations of the magnetic Sachs form factor is the same as the one for Ga and Gp. 



3. Results 



We perform the calculation at three pion masses with a fixed lattice spacing. The source and 
sink time slices are taken at = and tf = 12 and a check with tf — ti = 14 on 96 configurations 
at the highest mass gave compatible results, as shown in the right panel of Fig. The remaining 
parameters of the calculation together with the most important results are summarized in Table. ||. 



mjc[GeV] 


number of confs 


niN [GeV] 








gA 


0.4470(12) 


346 


1.287(13) 


0.489(25) 


0.558(61) 


2.83(21) 


1.171(41) 


0.3903(9) 


184 


1.245(9) 


0.527(34) 


0.607(79) 


2.90(36) 


1.096(47) 


0.3131(16) 


419 


1.143(11) 


0.649(38) 


0.63(11) 


2.85(46) 


1.23(10) 



Table 1: The form factors are extracted using a lattice of size 24^ x 48 at j3 = 3.9 with a lattice spacing of 
a = 0.0889(12). This table summarizes the main results at the three different pion masses. 



3.1 Electric and magnetic Sachs form factors 



The dependence of Ge and Gm on the euclidean momentum transfer squared is shown in 
Fig. m A dipole form 



Ge{Q^ 



1 + 



_Gm(0) 
1 + ^ 



(3.1) 



describes our data very well and the solid lines in Fig. |2| are the corresponding least squares fits. 
From the slope at zero momentum transfer, an electric and magnetic radius can be defined 



{''em) 



dGEM^Q^ 



12 



22=0 



m 



(3.2) 



£,M 



These can easily be translated into the more common Dirac and Pauli radii 



in) = -T 



12 3Gm(0)-3 



2rrfi 



Gm{^) ~ 1 V fi 



2Gm(0) 



+ 



2m?- 



(3.3) 
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, |J = 0.0085, m = 447 MeV 
y \i = 0.0064, = 390 MeV 
, H = 0.004, = 313 MeV 
experiment 
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Figure 2: Q^-dependence of Ge (left) and Gm (right). The dashed curve is a dipole fit to experimental data. 
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Figure 3: Comparison of results on (rj), {rf) and jj.iv between A^f = 2 twisted mass and Wilson fermions. 



Using the value of the lattice spacing determined from an analysis of the nucleon mass within 
the twisted mass QCD framework we convert our lattice quantities to physical units and give the 
results in Table ||. The magnetic form factor at zero momentum transfer is directly related to the 
magnetic moment, Gm{0):^ = fXiy. In Table |T] we give the iso vector magnetic moment converted 
to nuclear magnetons using the physical nucleon mass. We show the dependence of {r\), {r\) and 
jiiv on the pion mass in Fig. |3[ For comparison we include the corresponding results of Ref. [^], 
which were obtained with two degenerate dynamical Wilson fermions in a similar setup. As can 
be seen, the results in the two formulations are in good agreement. With the isovector current, 
the corresponding physical observable is the magnetic moment of the proton minus the one of the 
neutron. This is the experimental point shown in the right panel of Fig. |[ As can be seen, the 
experimental values in all cases are higher than lattice results. A chiral extrapolation of our data 
to the physical point will be carried out once we obtain results at an additional lighter pion mass. 
This additional input is needed in order to obtain reliable results at the physical point. 



3.2 Axial form factors 

The axial form factors are shown in Fig. ^with Gp showing the larger statistical errors. We get 
a good signal for the axial charge gA = ZaGa(O), with the renormalization constant Za = 0.76(1) 
that has been computed in [|l^]. We give our values in Table. ^. 
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Figure 4: 2^ -dependence of Ga (left) and Gp (right). The dashed green curve (left) is a dipole fit to 
experimental data on Ga yielding niA = 1.1 Ge V. In the right figure it derives from the fit to experimental 
results on Ga together with the assumption of pion pole dominance. 

4. Conclusions 

The isovector electromagnetic and axial nucleon form factors are evaluated using two dynam- 
ical degenerate twisted mass fermions for pion masses down to about 300 MeV. The results are in 
agreement with previous lattice studies [^. Like in Ref. we find better agreement with experi- 
ment in the case of the electric form factor than for the magnetic one. Although our results on the 
axial charge are in good agreement with experiment, we find a weaker momentum dependence for 
Ga- At low Q^, Gp shows a steep increase as expected from pion pole dominance and approaches 
the theoretical prediction as the pion mass decreases. In the future we plan to analyze configura- 
tions at a lighter pion mass and larger volume. This will enable us to reliably carry out a chiral 
extrapolation of the magnetic moment, radii and axial charge to the physical point. 
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